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We present data analyses of the thermal conductivity and specific heat of electron-doped 
Pri.85Ceo.i5Cu04_ !/ . The zero-field thermal conductivity and specific heat of this optimally electron- 
doped system can be only explained by a nodeless gap symmetry. The magnetic-field dependence of 
the electronic specific heat in the vortex state is in quantitative agreement with an s-wave theory. 
Our quantitative data analyses provide bulk evidence for a nodeless gap symmetry in optimally 
electron-doped cuprates. 



An unambiguous determination of the gap symmetry 
in the bulk of cuprates is crucial to the understand- 
ing of the pairing mechanism of high-temperature su- 
perconductivity. Three major gap symmetry contenders 
have been s-wave, d-wave, and extended s-wave (s + g 
wave). Both d-wave and extended s-wave have line 
nodes and change signs when a node is crossed. The 
majority of measurements probing low-energy excita- 
tions in the superconducting state (e.g., magnetic pen- 
etration depth, thermal conductivity, and specific heat) 
have pointed towards the existence of line nodes in 
the gap function of hole-doped cuprates Quali- 
tatively, these experiments are consistent with both d- 
wave and extended s-wave gap symmetries. In contrast 
to hole-doped cuprates, the earlier Raman scattering 
5] and magnetic penetration depth data [6) suggested 
a nodeless s-wave gap symmetry in optimally electron- 
doped (n-type) cuprates. Later on, phase and surface- 
sensitive experiments 0, H| provide evidence for pure d- 
wave symmetry in optimally doped and overdoped n-type 
cuprates. Surface-sensitive angle-resolved photoemission 
spectroscopy (ARPES) d [ijj also supports d-wave gap 
symmetry for optimally electron-doped cuprates. On the 
other hand, nearly bulk-sensitive point-contact tunneling 
spectra [llj show no zero-bias conductance peak (ZBCP) 
in optimally doped and overdoped samples 13l4l6j while 
the ZBCP is seen in the tunneling spectra of a deeply un- 
derdoped n-type cuprate with T c = 13 K (Ref. [l3lll3|V 
These tunneling experiments indicate a possible crossover 
from d wave in deeply underdoped materials to node- 
less s wave in optimally doped and overdoped materials. 
Very recently, the bulk-sensitive Raman scattering data 
of Ndi.ssCeo.isCuO^j, can be quantitatively explained 
by an anisotropic s-wave gap with a large minimum gap 
of about 3 meV (Ref. 17]), in quantitative agreement 
with the earlier penetration depth data d, [ljj. Fur- 
ther, the bulk- and phase-sensitive data of impurity pair- 
breaking effect in optimally-doped Pri.855Ceo.i45Cu04 



unambiguously rule out any d-wave gap symmetry 17 1. 
Therefore, the surface-sensitive experiments probe a d- 
wave gap at any doping level while the bulk-sensitive 
experiments see an s-wave gap in optimally doped and 



overdoped n-type cuprates. This apparent contradiction 
can be naturally resolved if the surfaces are deeply under- 
doped so that the gap s ym metry is d-wave. Experiments 
on hole-doped cuprates |l9l . [2fj| indeed show that surfaces 
and interfaces are significantly underdoped. 

Here we provide quantitative analyses of the bulk- 
sensitive thermal conductivity and specific heat data of 



optimally doped n-type Pri.ssCeo.isCuO 



4-v 



We find 



that the zero-field thermal conductivity and specific heat 
of this optimally electron-doped system can be only ex- 
plained by a nodeless gap symmetry. The magnetic-field 
dependence of the electronic specific heat in the vortex 
state is in quantitative agreement with an s-wave theory. 
The thermal properties of the optimally electron-doped 
cuprates consistently provide bulk evidence for a nodeless 
s-wave gap symmetry. 

The low-temperature electronic thermal conductivity 
K e i and specific heat C e i can be used to distinguish be- 
tween a nodeless gap function and a gap function with 
line nodes. Specific heat experiments are insensitive to 
the phase of the gap, but can provide bulk information 
on the behavior of the density of states N(E) near the 
Fermi level Ep. The quantity C e i/T is proportional to 
N(E) averaged over an interval ksT around Ep- For a 
d-wave or extended s-wave symmetry, the gap vanishes 
on lines of nodes on the Fermi surface, therefore N(E) 
averaged over all directions in the reciprocal space be- 
haves as N(E) oc |S| . This implies C ei /T oc T for T « 
T c in zero magnetic field. In contrast, for a fully gapped 
superconductor, C e i/T ~ under the same conditions. 
Similarly, the quantity n e \ cx T at T << T c in zero mag- 
netic field for d-wave or extended s-wave gap symmetry 
while for a fully gapped superconductor, k e i/T ~ under 
the same conditions. 

More specifically, the quantity K e i /T in zero magnetic- 
field is directly related to the Fermi velocity vf and mo- 
mentum hkp in the nodal directions, and to the slope 
S = dA(0)/d8 at nodes, where A(0) is a gap function and 
9 is the angle measured from the Cu-0 bonding direction. 
The former two quantities can be obtained from angle- 
resolved photoemission spectroscopic data while the lat- 
ter one can be readily calculated from the gap function. 
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The residual thermal conduction is due to a fluid of zero- 
energy quasiparticles induced by the pair-breaking effect 
of impurity scattering near the nodes in the gap. Calcu- 
lations for the heat transport by nodal quasiparticles in 
two dimensions give a general expression [J] 
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where n/d is the stacking density of Cu02 planes, vi = 
S/(Tikp). One can readily show that S — 2Am for a 
simple d-wave gap function: A(6) = Am cos 29 (where 
Am is the gap along the Cu-0 bonding directions). Sim- 
ilarly, the low-temperature electronic specific heat in zero 
magnetic-field is proportional to T 2 , i.e., C e i — aT 2 , 
where a is a constant. 

Another way to test the gap symmetry is to study the 
electronic specific heat in the vortex state. According to 
a simplified argument [21], the energy of carriers circu- 
lating around a vortex is shifted by the Doppler effect. 
This shift has a dramatic effect on the density of states 
when the gap is small, so that the essential contribution 
comes from the vicinity of the nodes. In the low temper- 
ature limit, the density of states at Ep becomes propor- 
tional to \fB for one vortex so that C e i/T oc \/H. In 
an isotropic s-wave superconductor, no significant con- 
tribution should be expected at low temperatures from 
such a mechanism, but localized states in vortex cores 
contribute a term which is proportional to H/H c %, that 
is, C e i/T oc H. However, the above simple picture should 
be valid only near the lower critical field H c \ (Ref. 22]). 
By taking account of the core excitations from the bound 
states in a vortex core and the quasiparticle transfer be- 
tween vortices, Ichioka et al. [22j showed C e i(T,H)/T 
for an isotropic s-wave superconductor is proportional to 
H o.67 for H < Hc2 an ,i T = while C e i(T,H)/T for a 
e?-wave superconductor is approximately proportional to 
only for H < Q.lH c2 and T = 0. 

Figure la shows k/T as a function of T for an in- 
sulating parent compound Nd2Cu04. The data are from 
Ref. [23] . Low-temperature thermal conductivity data of 
several superconductors have been fitted by an emperical 
formula: k/T = A + BT P , where A and B are constants, 
and the exponent p < 2 (Ref. [24|). The constant A is 
equal to K e i/T in a e?-wave superconductor and A = 
for a gapped s-wave superconductor. For the insulating 
material Nd2Cu04 where K e i = 0, so the phonon contri- 
bution up to 0.5 K can be fitted by the emperical formula 
with A = and p = 1.49 (solid line). Fitting the limited 
data above 0.38 K (dashed linear line) gives A = 0.45 
mW/K 2 cm and p = 1.7. The inferred large value of A = 
0.45 mW/K 2 cm for this insulating compound is unphysi- 
cal, indicating that the fit to the limited data points leads 
to a wrong conclusion. 

Figure lb shows k/T as a function of T 1 " 7 for electron- 
doped Pri.85Ceo.i5Cu04_j, (T c = 21 K) in zero mag- 
netic field. The data are from Ref. [25J. Since the 
phonon contribution K p h/T of the related compound 
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FIG. 1: a) k/T as a function of T 1 ' 7 for an insulating 
parent compound Nd2Cu04. The data are from Ref. [23l |. 
The solid line is the fitted curve by an emperical formula: 
k/T = A + BT P with A = and p = 1.49. The fit to the data 
above 0.38 K (dashed line) gives A — 0.45 mW/K 2 cm and 
p = 1.7. b) k/T as a function of T 1 ' 7 for for electron-doped 
Pri.85Ceo.i5Cu04- u (T c = 21 K) at zero magnetic field. The 
data are from Ref. [25|]. The solid line is the fitted curve by 
the formula: k/T = A + BT P with A = and p = 1.49. The 
fit to the data above 0.33 K (dashed line) gives A = 0.50 
mW/K 2 cm and p = 1.7. 



Nd2Cu04 follows a single power law with the expo- 
nent of 1.49 in the whole temperature range from 0.07 
K to 0.5 K, one should fit the data by the equation: 
k/T = A + BT 1A9 . The best fit (solid line) leads to 
A = K e i/T = -0.036±0.015 mW/K 2 cm. The negligibly 
small value of K e i /T should provide compelling evidence 
for a nodeless gap symmetry, in agreement with other in- 
dependent bulk-sensitive experiments mentioned above. 
On the other hand, the fit to the data above 0.33 K 
(dashed linear line) yields A = 0.50 mW/K 2 cm and p 
— 1.7, which appears to be consistent with d-wave gap 
symmetry. Since the fit to the limited data points of the 
insulating Nd2Cu04 leads to an unphysical conclusion, 
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the similar fit to the limited data points of the super- 
conducting Pri.85Ceo.i5Cu04_j, may also give rise to a 
wrong conclusion. 
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FIG. 2: a) k/T as a function of T 1 ' 5 for a heavily overdoped 
Tl 2 Ba 2 Cu06+j / (T c = 15 K). The data are from Ref. [2(|. 
Since k p h/T is proportional to T 1,5 (see Fig. la above), we 
plot k/T as a function of T to clearly see if there is a low-T 
downturn in electronic thermal conductivity. Following the 
definition of the parameter r in Ref. [28|]. we calculate r — 7.6 
using the reported sample dimensions (the length — 0.3 mm, 
the cross-sectional area = 0.2 mm x 0.01 mm) 26], and using 
the electronic contact resistance = 0.2 Q. (Ref. 26]) and Kei/T 
= 1.39 mW/K 2 cm. b) k/T as a function of T 1 ' 5 for for an 
overdoped Lai. g3Sro.i7Cu04. The data are from Ref. [53]. We 
calculate r — 2.3 using the reported sample dimensions (the 
length — 0.97 mm, the cross-sectional area = 1.26 mmx0.212 
mm) [2a|, a nd using the electronic contact resistance = 0.01 
n (Ref. Hi) and K el /T = 0.5 mW/K 2 cm. 

There are some other possible explanations for the neg- 
ligibly small residual linear term within the context of 
d-wave gap symmetry. The first possibility is that d- 
wave nodal quasiparticles become completely localized 
below 0.5 K. However, the data in the magnetic field of 
13 T indicate [29| that the complete localization would 
occur well below 0.1 K since K e i/T is still substantial 



(0.28 mW/K 2 cm) at 0.1 K. This apparent contradiction 
makes this interpretation very unlikely. The second pos- 
sibility is that the measurement gradually ceases to de- 
tect the electronic heat current as electrons fall out of 
thermal equilibrium with the phonon bath when the tem- 
perature is below 0.5 K. Two facts lead us to rule out 
the electron-phonon decoupling [28] as the main mech- 
anism for the drop. First, identical data are obtained 
by sending heat directly through the electron system us- 
ing photons (with no thermal resistance) (2^. Secondly, 
if this model [28| were relevant, the downturn should be 
more pronouned in heavily overdoped Tl2Ba2Cu06+i, (T c 
— 15 K) than the overdoped Lai.g3Sro.i7Cu04. This is 
because the parameter r that sensitively controls the ex- 
tent of the low-T downturn [28| is calculated to be 7.6 
for heavily overdoped T^B^CuOg+y, which is a fac- 
tor of 3.3 larger than that (2.3) for Lai.83Sro.i7Cu04. 
However, there is no low-T downturn in Ti2Ba2Cu06+ y 
(Fig. 2a) while the significant drop starts below about 
0.26 K in Lai. 83 Sro.i7Cu04 (Fig. 2b). Therefore, it is 
hard to imagine that this model would be valid only for 
Lai.83Sro.i7Cu04 and Pri.ssCeo.isCuC^-y but not for 
Tl 2 Ba 2 Cu0 6 + v . 

In Figure 3a, we show the temperature dependence of 
specific heat for a Pri.gsCeo.isCuC^-y crystal (T c = 23 
K), at different mag netic fields (H || c). The data are 
taken from Ref. [30(. If the optimally doped cuprate is 
a gapped s-wave superconductor, C(T, H)/T should be 
expressed as C(T, H)/T = 7 (0, H) + f3T 2 , where the first 
term is the field-dependent electronic specific heat and 
the second term is the lattice specific heat due to acoustic 
phonons. It is clear that the specific heat data at different 
magnetic fields can be well fitted by the above expression. 
The average of the fitted j3 values is 0.25 mJ/mol K 4 , 
leading to a Debye temperature of 384 K. A finite value 
of 7(0, 0) is also found in hole-doped cuprates. and there 
are several proposals to explain its origin |3ll - l33j . 

On the other hand, if electron-doped cuprates are d- 
wave superconductors with line nodes, C(T,0)/T should 
be expressed as C(T, 0)/T = 7 (0, 0) + aT + /3T 2 , where 
the term aT arises from nodal quasiparticle excitation 
(see discussion above). In Fig. 3b, we fit the zero-field 
data by C(T, 0)/T = 7 (0, 0) + aT + f3T 2 . Although the 
three-parameter fit improves the quality of fitting slightly 
compared with the above two-parameter fit, the fitted 
parameter a is negative (— 0.25±0.05 mJ/mol K 3 ) and 
unphysical. This suggests that the value of a is negligibly 
small, in disagreement with pure d-wave gap symmetry 
with line nodes. 

In order to obtain the zero-temperature 7(0, H) more 
precisely, we fit the data of Fig. 3a using C(T, H)/T = 
7(0, H) + I3T 2 with a fixed /3 = 0.25 mJ/mol K 4 . Fig. 4 
shows [7(0, if)— 7(0, 0)1 1 in as a function of H/H r2 . Here 
H C 2 is taken to be 50 kOe (see Fig. 2a of Ref. [30() and 7at 
= [7(0, 56kOe)-7(0,0)] = 4.1 mJ/mol K 2 . It is striking 
that the field dependence of the electronic specific heat 
for this electron-doped cuprate is in quantitative agree- 
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FIG. 3: a) Temperature dependence of specific heat for a 
Pri.85Ceo.i5Cu04_ !/ crystal (T c = 23 K), at different magnetic 
fields (H || c). The data are taken from Ref. f3(|. The solid 
lines are linearly fitted curves, b) Temperature dependence 
of specific heat at zero magnetic field. The data are fitted by 
C(T,0)/T = 7(0,0) + aT + f3T 2 . The fitted parameter a is 
negative (— 0.25±0.05 mj/mol K 3 ) and unphysical, suggesting 
that the value of a is negligibly small, in agreement with a 
nodeless gap symmetry. 



merit with a theory based on isotropic s-wave gap sym- 
metry 22J. This excellent agreement clearly indicates 



that the optimally electron-doped cuprate is a gapped 
superconductor . 

Fitting the zero-temperature data and theoretically 
calculated values respectively in Fig. 4 by a power 
law yields an exponent of 0.64±0.04 for the data and 
0.663±0.008 for the theory. At finite temperatures, 
the exponent should increase with increasing temper- 
ature, as demonstrated in an s-wave superconductor 
NbSe2 (Ref. [34j). In order to further verify a node- 
less gap symmetry in the optimally electron-doped sys- 
tem, Fig. 5 compares [C(T,H) - C(T,0)]/j N T for a 
P r i.85Ceo.i5Cu04_ y crystal (T c = 22 K) with that for 
V3SL The data of Pri.ssCeo.isCuOzi-,, are reproduced 
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FIG. 4: Magnetic field dependence of [7(0, H) - 7(0, 0)]/7jv 
for the Pri.ssCeo.isCuCU-j, crystal (solid circles) along 
with the theoretical prediction (solid diamonds) based on 
isotropic s-wave gap symmetry [22| . Here H C 2 and 7jv for 
Pri.85Ceo.i5Cu04_ H are taken to be 50 kOe and 4.1 mj/mol 
K 2 , respectively. 
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FIG. 5: Magnetic field dependence of [C(T,H) - 
C(T,0)]/7jvT for a Pn.gsCeo.isCM^-y crystal with T c = 22 
K at T = 1.8 K (open circles) and for V 3 Si at T = 3.5 K 
(solid circles). Fitting the two sets of data respectively by a 
power law gives the same power of 0.78. 



from Ref. [35| and the data of VaSi are from Ref. [36| . 
The values of 7jv and H C 2 for the PCCO compound are 
taken to be 3.7 mj/mol K 2 and 43 kOe (see Fig. 2 of 
Ref. [1^1), respectively. The value of H c2 for the VaSi 
superconductor is found to be 250 kOe (Ref. (36|). Using 
the specific-heat jump at T c (i.e., AC/T C ~ 90 mj/mol 
K 2 
V 3 Si 



3611) and the relation AC/7atT c = 2.0 suitable for 
371 ]. we calculate 7Ar to be 45 mj/mol K 2 . It is ev- 



ident that the two sets of data almost coincide with each 
other. Fitting the two sets of data respectively by a power 
law produces the same exponent of 0.78, which is also 
very close to that (0.75) for NbSe 2 at 2.3 K (Ref. (H). 
Therefore, electron-doped cuprates have a very similar 
field dependence of the electronic specific heat as the con- 
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vcntional s-wave superconductors. 

In summary, independent bulk-sensitive specific heat 
and thermal conductivity experiments show a nodeless 
gap symmetry in the bulk of optimally electron-doped 
cuprates. A T 2 dependence of the penetration depth at 
low temperatures observed in some samples |38[ can be 
well explained [18| in terms of nodeless s-wave gap along 
with an extrinic effect due to current-induced nucleation 
of vortex-antivortex pairs at defects. Although the data 
of spin-lattice relaxation rate and Knight shift in an un- 
derdoped Pro.9iLaCeQj)gCu04_j, are consistent with d- 
wave gap symmetry [39j|). we will show that the data can 
be better explained by an anisotropic s-wave gap with 
a minimum gap size being consistent with independent 
magnetic penetration depth data. Further, the d-wave 
gap symmetry in Pro.giLaCeo.o9Cu04_j, is incompatible 
with the observed large residual resistivity (92 /iSlcm) 
[39I ] . which would suppress T c to zero [13] ■ Therefore, 
these bulk-sensitive experiments consistently point to- 
wards a nodeless s-wave gap symmetry, in agreement 
with predominantly phonon-mediated pairing mechanism 

us 53. 
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